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1. Introduction 

The quantitative understanding of the phase diagram of nuclear matter is a very important open 
problem which defies both analytical and numerical treatments [1], On the analytical side, sketches 
of the QCD phase diagram can be drawn from mean field arguments and the analysis of simple 
truncated models, which are mostly conjectural. On the numerical side, the notorious sign problem 
prevents the use of standard Monte Carlo methods to simulate lattice QCD at finite baryon density 
from first principles. 

Traditionally, the fermion fields in the partition function of lattice QCD are integrated analyt¬ 
ically a priori, which reduces the system to an ensemble of color-charged gauge fields weighted 
by the determinant of the lattice Dirac operator. The sign problem occurs for finite values of the 
baryon chemical potential, for which such a determinant is, in general, complex-valued. The direct 
Monte Carlo sampling of this ensemble is therefore impossible. Several numerical methods have 
been developed to deal with this difficulty, but none has succeeded yet in providing a satisfactory 
description of the QCD phase diagram. 

The common denominator to all such methods is the a priori integration over the lattice 
fermions, and the subsequent sampling over degrees of freedom carrying color-charge. However, 
QCD is known to be confining, in the sense that its asymptotic eigenstates are color-neutral. Hence, 
it would be more “natural” for the Monte Carlo integration to sample over color-neutral degrees of 
freedom - at least in the confining regime - rather than sampling over local gauge-covariant states, 
whose quantum fluctuations and local gauge ambiguity propagate into uncontrollably large fluctu¬ 
ations of the phase of the Dirac determinant. More precisely, if |!//,} denotes the putative complete 
basis of eigenstates of the Euclidean QCD Hamiltonian //, 



( 1 . 1 ) 


then all the matrix elements {\jfi\e~^\\jfj) are non-negative by definition, hence the sign problem 
disappears in such a basis. In short, the severity, and even the existence, of a sign problem depends 
on the particular representation of the partition function, Z. 

Even though it is unfeasible to represent the partition function of lattice QCD in terms of its 
physical eigenstates, it is reasonable to think that the sign problem might be significantly reduced 
when Z is represented in terms of color-neutral states which are not necessarily eigenstates. Such 
a representation can be obtained by integrating the lattice gauge fields a priori. Let us consider the 
partition function of SU (N) or U (N) lattice gauge theory with a single flavor of staggered fermions, 

2= j\dx^x\ ^ Tp(t — yyReTr(t/p))^£ v ^ il.\iiT'r(XxUxnXx-t-fi—Xx+fL^lxiiX-x)^ am Y.xXxXx (| 2) 

where U xp are link variables, U p is the path-ordered product of link variables around a plaquette, 
XxiXx are staggered fermions, am is their hare mass in lattice units, /3 is the bare lattice coupling, 
and T/ xp = ±1 are staggered phases. The direct analytic integration of the gauge fields is only 
feasible at /3 = 0: the plaquette terms vanish, and the fermionic weight factorizes into a product of 
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Figure 1: Typical MDP configuration of Nf = 1 lattice QCD on a 6 x 6 lattice, at /3 = 0. 
solvable one-link integrals [2], 

Z = / [dxdx\e 2am ^^\\ f (1.3) 

J X,/A J 

= I [dxd X ]e^U (jt + (1.4) 

The resulting partition function is a Grassmann integral whose integrand contains color-neutral 
terms only, namely mesons M x = XxXx and baryons B x = jfi£i v ..i N Xx '' 'X[ N ■ Color-neutral states 
made of such building blocks are not exact eigenstates of lattice QCD, but they resemble the ex¬ 
pected asymptotic states of the continuum theory. 


2. Dimer representation of lattice QCD 


The lattice fermions in Eq. (1.4) still need to be integrated out. Solving the Grassmann integrals 
results in a combinatorial partition function [2], 


Z = 


I 

{n,k,(} 


n^w-n 

X,g 


x n x- 


{N-k xfl )\ a(£) 
N\k xll \ AHM 


(2.1) 


whose degrees of freedom are monomers, dimers and self-avoiding polygons which saturate the 
spacetime lattice. More precisely: 


• n x G {().... ,/V} represents the number of monomers occupying the lattice site x. It is in¬ 
terpreted as the condensate contribution (M x ) nx to the fermionic state. The weight of such 
monomers depends on the bare quark mass am and is positive-definite. 

• k x n G {0,... ,N} represents the number of dimers occupying the lattice bond (jc,/t). It is in¬ 
terpreted as the meson hopping contribution (M x M x+ h) kx > 1 to the fermionic state. The weight 
of such dimers is also positive-definite. 

• l is an oriented self-avoiding polygon. It is interpreted as the contribution to the fermionic 
state of a sequence of baryon hoppings B x B x+ fi which form a closed loop. The weight of a 
baryonic loop t depends on its perimeter |f| and on its sign a(£) = ±1, which is a specific 
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Figure 2: Graphical representation of the Q variables (left) and of the R variables (right). 

function of the topology and shape of l. The sign problem occurs due to the existence of 
baryonic loops £ for which a(£) = — 1. 

The Grassmann integration constrains the monomer-dimer-polygon (MDP) configurations, 
only allowing those which saturate each site with exactly N quarks and N anti-quarks. Each bary¬ 
onic loop must therefore be a self-avoiding polygon which saturates every site it crosses, while 
the number of monomers or dimers attached to the remaining sites must be exactly N = n x + 
£ ±jU k xfl , fix (see Fig. 1 for example). 

Constrained combinatorial systems such as Eq. (2.1) can be simulated efficiently with variants 
of the worm algorithm [3]. Recent simulations of the MDP representation of lattice QCD (N = 3) 
using such algorithms [4] show that the sign problem becomes very mild, with an improvement by 
a factor of 0( 10 4 ) in accessible volumes over the traditional approach. The sign problem is also 
shown to remain mild when 0(j8) corrections are included. Consequently, the phase diagram of 
lattice QCD in the strong coupling limit and its 0(/3) corrections have been mapped successfully 
[4]. 


3. Decoupling the plaquette 


Higher-order corrections to the strong coupling limit of lattice QCD in the MDP representa¬ 
tion involve combinatorially large numbers of cumbersome group integrals. In order to approach 
the regime of continuum physics, it would be desirable to use a MDP representation of lattice QCD 
which is not of the form of a strong coupling expansion. However, it is not possible to directly 
integrate out the link variables in the presence of Wilson action terms, which couple the four link 
variables around each plaquette. We circumvent this analytical difficulty by coupling the link vari¬ 
ables to auxiliary fields in such a way that that links decouple among themselves [5]. We illustrate 
this in the case of the pure gauge theory, 

Z = |[Jf/] e -^ P ( 1 -^ ReTr M (3.1) 


We refer to Eq. (3.1) as the 4-link partition function, given that it only contains terms with four link 
variables coupled around a plaquette. 

First, we introduce a set of free auxiliary variables Q XJU v £ C NxN living on “diagonal” links 
(see Fig. 2), i.e. with Gaussian distribution dQ x/xv dQ x ^ v exp(—^Tr(Q xluv Q x ^ v )). A Hubbard- 
Stratonovich (HS) transformation of the form [5, 6]: 


Qxgv ^ 



+ UjcgUx+fi^ 


T U xv U x ^y fj ) 


(3.2) 
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Figure 3: Graphical representation of the HS transformation of the Q variables (top), which induces the 
replacement of the Wilson plaquette terms in the 4-link partition function by j-plaquette terms (bottom). 

splits each Wilson plaquette term into two if-plaquette terms, each of which couples two link vari¬ 
ables U and one auxiliary variable Q (see Fig. 3). The partition function Eq. (3.1) then becomes: 


Z« / [dQdQ'je 


? Tr (ee+ 


3 N vEvaG-u.v)) 


(3.3) 


which we refer to as the 2-link partition function, given that it only contains terms with two link 
variables coupled around a 4-plaquette. 

Second, we introduce yet another set of free auxiliary variables R x ^ v e C' Vx,v (i.e. with Gaus¬ 
sian distribution) living on “folded” links (see Fig. 2). A HS transformation of the form [5]: 


R 


XflV 


^ \l Xji R xilV + Qxtl V U l + fiy + U xg) 


N 


(3.4) 


splits each ^-plaquette term of Eq. (3.3) into two ^-plaquette terms, each of which couples a single 
link variable U to ambient auxiliary variables Q and R (see Fig. 4). The partition function Eq. (3.3) 
then becomes: 


Z« [[dRdR^ldQdQ^e-^™ 1 ^-^ 00 ^ 


J] f dUe* ReT < J » 




where J xfl is a function of the auxiliary variables attached to the link (. x,fi ), 


Jxg — ^ y^x-v,vgQx-v,vg,~r^xgv 

vfig 


+Rx 


(3.5) 


(3.6) 


We refer to Eq. (3.5) as the 1-link partition function, given that it only contains ^-plaquette terms 
consisting of a single link variable coupled to its ambient auxiliary fields. 

The integrand of Eq. (3.5) factorizes into a product of one-link integrals, which can be solved 
analytically [7]. For example, the 0-link partition function of compact U (1) lattice gauge theory is 
given by: 


j [dRdR*][dQdQt]e~% 


rr\- 3 4q& 


iimm 

x,g 


(3.7) 


where I v (z) are modified Bessel functions of the first kind. We refer to Eq. (3.7) as a 0-link partition 
function, given that all link variables have been integrated out. 
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Figure 4: Graphical representation of the HS transformation of the R variables (top), which induces the 
replacement of the ^-plaquette terms in the 2-link partition function by ^-plaquette terms (bottom). 


For the various //-link representations with n > 0, loop observables can be computed ordinarily 
as path-ordered products of link variables around the loop £:U(£) = Ylirt U\■ In the 0-link represen¬ 
tation, however, the link variables are also integrated out in the definition of the loop observables. 
It is straightforward to show that the 0-link loop observables are path-ordered products of effective 
links U xfl :U(£) = Ui e t U/, so that 

(±T*(U(£))^ = (±Tr(U(£))) (3.8) 


The effective links U x ^ are defined as the first moments of the corresponding one-link integrals: 


Ug = (U ij ) Jx)1 = ^f dUU ij e* ReT < J *^ 

In compact (7(1) lattice gauge theory, for example, the effective link is simply given by: 

h (P\Jxg\) J m 


U m = 


UP M) M 


(3.9) 


(3.10) 


In Monte Carlo simulations of the //-link partition functions, for // > 0, the U, Q,R variables 
are updated as follows: 

1. Update Qxnv via Gaussian heatbath with variance j and mean U x ^U x+ f l)V + U xv U x+ ^^\ 

2. Update R x gv via Gaussian heatbath with variance j and mean Q X ^ V U^ + ^ v + ; 

3. Update U x ^ via (pseudo)heatbath with respect to the Boltzmann weight e ^ ReTl ^ Ux ^ 


The same algorithm can be used to simulate the 0-link partition function, in which case the link 
variables U play the role of auxiliary variables, while the Q. R play the role of dynamic variables; 
and so the U variables must then be updated prior to the Q. R variables. 

Numerical simulations confirm that the expectation value of loop operators (e.g. the plaquette 
operator) is reproduced exactly in all //-link representations, in the entire range of (5. However, the 
more auxiliary fields are present, the larger are the errors associated with lattice measurements, and 
the larger are the corresponding autocorrelations. In short, the //-link representations with // < 4 arc 
not adequate for simulations of the pure lattice gauge theory. 
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4. Compact lattice QED 


Given that the staggered fermion action is linear in the link variables, the 1-link representation 
Eq. (3.5) of the pure gauge theory can be used to extend the MDP representation of the strong 
coupling limit Eq. (2.1) to arbitrary values of /3, by shifting / vU before performing the group inte¬ 
gration: 

2N Jxv ^ 2N Jxv + Vxv^^XxXx+v (4-1) 

2 N J tv ^ 2N J x y ~ (4 ‘ 2) 

Here pi denotes the quark chemical potential. 

We illustrate this extension with Nf = 1 compact lattice QED on a lattice torus [8], for which 

the electron chemical potential vanishes due to Gauss’s law. Its 0-link partition function at finite f5 

is given by: 


2 2/ [dRdR'}[dQdQ f ] e-2 RR \- 3 i Q(I ]“[ / dU 
J x,n J 


j [, dRdR r ] \dQdQ 






x,g 


g 2 amZxXxXx x 


x (1 + M x M x+(l + ti^U^XxXx+h ~ flxgUxgXx+piXx) 


(4.3) 


In addition to the “meson” hoppings also present in the strong coupling limit, Eq. (4.3) contains 
U (l)-invariant electron hoppings of the form U^XxXx+p.^ which vanish at (5 = 0. After Grassmann 
integration, we obtain the MDP representation of compact lattice QED at finite f5: 

Z oc j[dRdR^}[dQdQ^]e~ We-ZottllWM) £ l\(2am) n *o(£)U\£) (4.4) 

4 x,fl {n,k,e} x 

°< £ Y\^c{£){U\£)) (4.5) 

{n,k,£} x 

In the strong coupling limit, Eq. (4.5) reduces to Z = Tlx (2ara)”*, which sums over monomer- 
dimer coverings of the spacetime lattice. At finite /3, in addition to monomers and dimers, there are 
also contributions from self-avoiding electron loops £, whose weight is given by the pure-gauge 
expectation value of the Wilson loop around £, multiplied by an overall sign (see Fig. 5). Such a 
representation corresponds to an exact loop expansion of the determinant of the Dirac operator. 

This representation has a sign problem due to the sign a(£) of the electron loops and the sign 
of Re ((/(£)). In the strong coupling limit, the sign problem vanishes and is therefore expected to 
be mild at small values of /3, where the typical configuration is a dilute gas of small electron loops 
in a monomer-dimer bath. For increasing /3, however, the entropy of electron loops dominates, 
and consequently the system is expected to enter a regime in which the sign problem is severe. 
Nonetheless, the sign problem is expected to be manageable in some window of small /3. 


5. Summary 

We showed that it is possible to fully decouple the link variables around a plaquette in the 
Wilson gauge action, at the cost of coupling them to a new set of free auxiliary fields. Consequently, 
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Figure 5: Typical MDP configuration of Nf = 1 compact lattice QED on a 4 x 4 lattice, at finite /j: the red 
“diagonal” links and the blue “folded” links represent the auxiliary Q and R variables, respectively; dots 
represent monomers, solid bonds represent dimers, and dashed lines represent electron loops. 

the gauge fields can be integrated out exactly for any value of /3, thus resulting in a representation 
of the lattice partition function without link variables. In the 0-link representation, however, errors 
are larger and autocorrelation times are longer than in the original 4-link representation, which 
renders it unattractive for simulations of the pure gauge theory. 

Nonetheless, it allows us to integrate out the link variables exactly also in the presence of 
an arbitrary number of (staggered) fermion flavors. We illustrated this with the simple example of 
compact lattice QED on a torus, for which we obtained an exact MDP representation of its partition 
function at arbitrary /3, as a sum over electron loop configurations. Numerical simulations using 
worm algorithms are in progress. 

We also plan to construct MDP representations of lattice gauge theories with non-trivial de¬ 
pendence on the chemical potential, e.g. Nf > 1 lattice QED, and non-Abelian gauge theories 
including lattice QCD. The corresponding partition functions are obviously more complicated, due 
to the proliferation of gauge-invariant terms that contribute to their 0-link partition functions. But 
the recent successful simulations of the strong coupling regime of lattice QCD at finite density in 
the MDP representation [4] suggests that a larger window of lattice couplings might be accessible 
to Monte Carlo simulations using the methods presented in this talk. 
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